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SPATIAL EXTREMES: MODELS FOR THE STATIONARY CASE 

By Laurens de Haan 1 and Teresa T. Pereira 2 

Erasmus University and University of Lisbon 

The aim of this paper is to provide models for spatial extremes 
in the case of stationarity. The spatial dependence at extreme levels 
of a stationary process is modeled using an extension of the theory of 
max-stable processes of de Haan and Pickands [Probab. Theory Re- 
lated Fields 72 (1986) 477-492]. We propose three one- dimensional 
and three two-dimensional models. These models depend on just one 
parameter or a few parameters that measure the strength of tail de- 
pendence as a function of the distance between locations. We also 
propose two estimators for this parameter and prove consistency un- 
der domain of attraction conditions and asymptotic normality under 
appropriate extra conditions. 

1. Introduction. The paper develops a framework as well as concrete 
models for statistics of spatial extremes which are sufficiently simple to be 
used in applications. Only the case of stationary processes is considered 
and the dependence structure will be represented by one parameter or a 
few parameters. Instead of developing more complicated models we aim at 
developing several simple models with somewhat different features. 

For simplicity of exposition and in order to stay close to the existing 
literature we shall start discussing processes which are defined on R rather 
than M?. After that we discuss stationary processes in R 2 . 

The setting is as follows. Consider independent replications of a stochastic 
process with continuous sample paths 

{X n (t)} te R, 
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n = 1, 2, Suppose that the process is in the domain of attraction of a max- 

stable process, that is, there are sequences of continuous functions a n > 
and b n such that as n — ► oo 

n f maxi<i< n Xi(t) -6 n (t) | if ; - / u 

i 1 - 1 ) 1 ttx f -n^wheR 

in C-space. Necessary and sufficient conditions have been given by de Haan 
and Lin [4]. The limit process {Z(t)} is a max-stable process. Without loss 
of generality we can assume that the marginal distributions of Z can be 
written as 

exp{-(l + 7 (t)x)- 1/7(t) } 

for all x with 1 + j(t)x > where the function 7 is continuous. For the time 
being we shall discuss the standardized process, called simple max-stable, 

{Z(t)} tm :={(l+ 7 (t)Z(t))^\ m , 

whose marginal distribution functions are all Frechet: exp(— 1/x), x > 0. 

{Z(t)} is a simple max-stable process. We assume that {Z(t)} is a sta- 
tionary process. The theory of de Haan and Pickands [5] applies. According 
to Theorem 6.1 of that paper the process is determined by a nonnegative 
L\ function and a group of linear L\ + isometries. However, since we aim at 
manageable models, we shall restrict ourselves to the subclass of stationary 
max-stable processes which is discussed on pages 490-491 of [5], the one of 
moving maximum processes. The process is defined as follows. 

Let (j) be a unimodal continuous probability density on the real line and 
let {Xj,Yj}j>i be the points of a homogeneous Poisson process on R x M + . 
The process is defined as a functional of the Poisson process as follows: 

Z(t) : = max ^ Xj - Q for t G R. 

3>l Yj 

It is easy to check directly that this process is stationary and simple max- 
stable. The almost sure continuity of the process follows from [1]. We think 
of t as a space parameter, not a time parameter. Note that for t±, £2, • ■ • , td £ 
M and X\,X2-, ■ ■ ■ , Xd > (cf. [5]), 

P{Z(t\) < x\, . . . , Z(tfi) < xA = exp{— / max — — ds\. 

I J-00 l<i<d Xi J 

However, this is not yet sufficiently simple for applications. We shall con- 
sider three specific examples depending on just one parameter. For <fi we 
choose the normal density 

(1.2) -^exp 
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the double exponential density 



(1.3) 



-exp{-/3|x|} 



and the i-density 



^T{u/2) 



{ 




(u+l)/2 



(1.4) 



1 + 



with v a positive integer 



where (3 is a positive constant. The constant /? measures the strength of tail 
dependence. We shall see that in all cases small values of (3 point at strong 
dependence and large values of (3 point at weak dependence. However, the 
dependence does not decrease at the same rate in all cases. 

As it should be for spatial models, the tail dependence between Z(0) and 
Z(t) decreases monotonically and continuously with \t\. In particular, when 
\t\ — > oo, the random variables Z(0) and Z(t) become independent. 

The same happens for fixed t when varying j3: as (3 I the process be- 
comes a.s. constant and as (3 — > oo Z(0) and Z(t) become independent (cf. 
propositions below). The dependence decreases monotonically as (3 increases. 

Next we extend the definition to processes on R 2 , that is, to random fields. 
It is readily seen that the theory of de Haan and Pickands [5] remains valid if 
the underlying Poisson process is based on M 2 x R + rather than R x R+. So 
we consider a unimodal (i.e., nonincreasing in each direction starting from 
the mode) continuous probability density <fi on R 2 . Let {Xj,Wj,Yj}j>i be 
the points of a homogeneous Poisson point process on R 2 x R + . The process 
defined by 



is easily seen to be stationary and simple max-stable. The a.s. continuity of 
the process follows from an extension of the arguments in [1]. 

The specific models we consider are analogous to the ones in the one- 
dimensional situation: 



Z{h,t 2 ): 



= max 



fijXj-tuWj-h) 



for (*i,t 2 ) 6l 2 



(1.5) 




(we call this the normal model), 



(1.6) <P(h,t 2 ) = ^e W {-P(\h\ + \t2\)} 



(we call this the exponential model) and 



(1.7) 




a > 1 
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(we call this the i-model), for /3 > 0. Finally we consider the general normal 
model 

(1.8) 4>{h,h) = ~ y= 2 exp(- . n 1 2 A P!tj - 2 P p 1 /3 2 tit2 

2n y / l-p 2 I 2(1- p 2 ) J 

where p is the correlation coefficient (— 1 < p < 1), for Pi, 02 > 0. 

The paper is organized as follows. The two-dimensional distributions of 
the process are derived for the mentioned models in Section 2. Those are suf- 
ficient for the estimation theory developed in Section 3. The two-dimensional 
marginal distribution for the normal model was derived earlier by Smith in 
an unpublished paper [14]. Higher-dimensional marginal distributions do not 
seem easy to calculate explicitly. 

In Section 3 we are mainly concerned with estimating the dependence 
parameter (3 on the basis of observations at finitely many locations from a 
stationary stochastic process in the domain of attraction of the max-stable 
process. In all the models there is a simple relation between (3 and a well- 
known dependence coefficient for two-dimensional extremes, 

A := limt -1 P{l - Fi(Xi) < t and 1 - F 2 (X 2 ) < *}, 
tjo 

where (X\,X2) has a distribution function F which is in the domain of 
attraction of some extreme-value distribution (Pi and F2 are the marginal 
distribution functions). 

The coefficient A 6 [0, 1] is related to the general framework of multidi- 
mensional extremes in the following way. If 

lim F n (a n x + b n ,c n y + d n ) = G(x,y), 

n— >oo 

where the two marginals of G are of the form exp{ — (1 + 7ix) + 1//7l j, i = 1, 2, 
then 

limt -1 P{l - FUXi) < tx or 1 - F 2 (X 2 ) < ty} 

V 71 72 ) 

= : L(x,y) 

and 

limt _1 P{l - Pi(Xi) < tx and 1 - F 2 (X 2 ) < ty} 
= x + y - L(x, y) =: R(x, y). 



Then A = 12(1,1). 
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D. Mason and X. Huang (see [9]) proved consistency and asymptotic nor- 
mality for a natural estimator R(x,y) of R(x,y). We use this for estimating 
A = R(l, 1). This result leads to a consistent asymptotically normal estima- 
tor of (3 based on observations taken at just two sites, £, and tj. In general 
observations are available at sites h,t2, ■ ■ ■ ,td (i.e., finitely many). One of 
our estimators of (3 is the average of /3-estimators based on the various pairs 
of sites. Consistency and asymptotic normality follow. 

The theory we develop can be used to solve some common problems in 
spatial extremes. One is an extrapolation problem: it consists of estimating 
the extremal behavior of a process X(t) at a site to where no observations 
are available based on repeated observations of the process at d different 
space points t±, tz, . . . , t^. Another problem refers to the extreme values of 
the (unobserved) aggregate process J s X(t) dt over a space region S assuming 
that the process has been discretely observed at a number of space points 
in S. Similarly one can look at the tail behavior of sup tgS X(t). 

In order to attack those problems, the estimation of (3 has to be com- 
plemented by estimation of local parameters: the extreme-value index, the 
scale and the location. The latter objects are not needed for the estimation 
of (3 but only for the application. 

The proposed models are quite simple examples of the representation in 
[5], which is valid for all stationary max-stable processes. It seems that the 
full model is not easily applicable: how does one estimate the initial non- 
negative L\ function and a general group of transformations? Moreover the 
representation is not unique. Instead we have tried to look at simpler mod- 
els which can be analyzed mathematically. We hope that providing several 
simple models with quite different features will widen the scope for applica- 
tions. Later we want to consider somewhat more general parametric groups 
of transformations. 

The validity of the model in applications can be checked with the following 
steps: estimate (3 as outlined in Section 3, estimate R(l, 1) for each pair of 
sites using the relation in Corollary 3.2 and check if this estimate of R(l, 1) 
is similar to the direct estimate using Proposition 3.2. We have not done 
this yet. 

2. Marginal distributions for the two models, and in M. and R 2 . We find 
the two-dimensional marginal distributions for all the models introduced in 
Section 1. Note that in this section, different from other sections, for sim- 
plicity and without loss of generality we consider the standardized process 
Z, not Z. First we consider the processes on the line, next the ones on M. 2 . 

Proposition 2.1. For t e M and the exponential model, for w\,W2> 0, 
-logP{Z(0)< Wl ,Z(t)<w 2 } 
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a roo r P -P\s\ p -p\s-t\ ) 

(2.1) =£/ max ,- ids 

— , for <W2 < e'P^wi, 

i i p -/3|t|/2 

(2.2) =<— + ■ : /or e"^V < w 2 <e^wi, 

W± W 2 y/WlW2 

— , for W2 > e^wi, 



(2.3) 



/ max< , >s{9)d0 

Jarctan(e-^l*l) I »1 »2 J 

1 111 f 1 e^l'l 

H — max< , — > H — max< 



where the spectral density s is given by 

P -P\t\/2 

(2.4) (fl)= - (singcos^)' 3 / 2 ; 

1 1 1 (w 2 \ 

= — + — + —x( — ), 

W± W2 W2 \WlJ 

and where the dependence function x ^ s given by 

!—s, < s < e~^l'l, 

-e~PW^ t e -/»l*l< a <e^l, s > 0. 

-1, s>e^l*l, 

Remark 2.1. Formula (2.1) follows directly from [5] and (2.3) reveals 
the spectral measure, which has a density on the interval (arctan(e~^'*'), 
arctan(e^'*')) and atoms of size at each of the two bound- 

ary points of that interval. The last characterization (2.4) is in the spirit of 
Sibuya [13] and Pickands [11]. 

Proof of Proposition 2.1. The integrand of (2.1) is f^p-^ if 

> £*!=L, that is, if 

inn tm 1 ' 



(2.5) | a _ t |_| s |<^iogf— 

P \W2 

Since the joint distributions of (Z(0),Z(t)) and (Z(0),Z(—t)) are the same 
we proceed as if t is positive. The left-hand side is t for s < 0, t — 2s for 
< s < t and —t for s > t. Hence if ^log(^) > t, then inequality (2.5) 

holds for all s and we get the first line of (2.2). Similarly if ilog(^) < -t, 
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we get the last line of (2.2). Next suppose —t < -^log(^) < t. In this case 
(2.5) becomes 



1 , 

t - 2s < - log 

P \W2/ 



that is, 



t 1 fW! 

2 2(3 \w 2 
Hence the integrand over this interval becomes 

■ ds = / ds 



2 Jt/2-l/(2/3)\og(w 1 /w 2 ) W2 U>2 2 J_oo W 2 

= 1 3_ e -/3t e) 8t/2 e -(l/2)log(wi/«j 2 ) 

W2 2lV2 

1 1 e-^l/ 2 



W2 2 y/WlW2 

This in combination with the integral stemming from the case e — < 
6 ^" gives the second line of (2.2). 

To check the equivalence of (2.2) and (2.3) for e~^*l < ^ < eP^ it suffices 
to see that the density of (2.2), after transformation to the polar coordinates 

^l+^i an d $ = arc t an ^f') is r~ 3 s(8) (cf. the construction of the 
spectral measure in [6]). For (wi,W2) outside this range just evaluate the 
integral in (2.3). □ 

Remark 2.2. The parameter P controls the dependence: if f3 — > oo, the 
spectral density s(9) goes to zero and the spectral measure is concentrated at 
the points {0, |}. This means that X(0) and X(t) are independent. If ft J. 0, 
the spectral measure concentrates on {j}- This means that X(0) = X(t) a.s. 

Proposition 2.2. For tsK and the normal model, for Wi,w 2 > 0, 
-logP{Z(0) <wt,Z(t) <w 2 } 

1 -{0\t\ 1 . Wx\ 



^ 2 ' 6 ^ t«r\ 2 1 p\t\^ Wl j 1 w 2 ^y 2 1 p\t\~ a w 2 j 

(2.7) =/ max , }s(9)d0, 

Jo I Wi U>2 J 

where the spectral density s is given by 
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+ G + ^ ln(tan 0) ) K^r - ^ ln(tan ' 



with (f(u) = $'(u); 



, N 1 1 1 

2.8 = — + — + — x 

Wl W2 W2 \W\ 



( w 2 
— 



where the dependence function x is given by 

*M = - s *(-if + 5 >o. 

Remark 2.3. Again the parameter 3 controls the dependence: if 8 —* 
oo, the variables Z(0) and Z{t) become independent; if 8 J, 0, we get Z(0) = 
Z(t) a.s. 

Remark 2.4. This distribution function has been obtained in a number 
of ways in the literature. Eddy [7] found it when studying convex hulls of 
samples. Hiisler and Reiss [10] obtained the distribution as the limit distri- 
bution of the componentwise maxima in a triangular array where the dis- 
tribution of the nth array is the two-dimensional normal distribution with 
correlation coefficient p(n) such that lim n ^ 00 (l — p(n))logn exists. A re- 
lated reference is [2]. In [14] Smith developed the distribution in the same 
way as it is done here. The distribution is mentioned in [3] and [12]. Falk, 
Hiisler and Reiss [8] obtained the distribution as the pointwise maximum 
of independent Brownian motions shifted by an amount corresponding to 
points of a Poisson point process. Another way of obtaining the distribution 
is 

-logP{Z(0)< Wl ,Z(t)<w 2 } 

= Emax{—, — exp(N8t-8 2 t 2 /2)\, teR, 

I W\ W2 ) 

where N is a standard normal random variable. 



Proof of Proposition 2.2. Clearly the distribution depends only on 
|t|. So we consider t > only: 

Q roc f e -/3 2 « 2 /2 e -/3 2 (u-t) 2 /2^ 

-logP{Z(0)<wi,Z(t)<w 2 } = -p= / max , \ du. 

V27TJ-0O I ^1 W 2 J 



Now e /3 '" 2/2 > e Phu t)2/2 if and only if Bu < + i In sa, hence 



logP{Z(0) <w 1 ,Z(t)<w 2 } 
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I l rf3t/2+{l/f3t)\nw 2 /w l 

e " 7 2 du 



W\ y/2TT J-oo 

1 1 



W 2 V27T J /3t/2+(l/l3t)\nw2/wi 



1 //W 1 , W 2 \ If fit 1 , Ttf 2 

— $ — + — In— + — -— H In — 

wi V 2 /ft wi 7 w 2 I V 2 /ft wi 



Hence the first part of the result. In order to obtain the second result note 
that 

' \ogP{Z(0)<w 1 ,Z(t)<w 2 }] 



(cf. the construction of the spectral measure in 



wi =r cos &,W2 =r sin t 



□ 



Proposition 2.3. For teR and £/ie t-model, for wi,w 2 > 0, 
-logP{Z(0) <W!,Z(t) <w 2 } 



(2.9) 



where 



1 



w 2 

— Pi,v(P,t,x) H (1 -p^C/MjS)), 

U>1 ^2 



-P< T Vtl < 



my 



0<w 2 <b^ 2 Wl , 



b 2 I wi < w 2 < wi, 



w 

(! -Pl,u{P^ X )) H P2,u(P,t,x), 

W\ w 2 



1 

U>1 



u>i = ^2 =: 



< ^2 < Ol^ Wi, 
. ,-(i/+l)/2 



p 2 Jp,t,x) 



1 + 




2zy 


1 + 


pH 2 


2u 


'{ 









P 2 t 2 



+ 



'v \ 

pt 

1-2 

/ftcc 

1 — 2 



1 + 



/? 2 t 2 



t 2 2 V 

(1-x) 2 

t 2 X f 
(1-X) 2 ~ ^ 
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T U) \ is a random variable with a Student t- distribution with v degrees of 
freedom and scale parameter 1 and x= (^) 2 ^ u+1 \' 



(2.10) 

where the dependence function \ is given by 



1 1 1 fW2 

W\ W2 W2 \WlJ 



X{s) 



0<s<b 



sp 1>v (P, t, S - 2 /("+D )-s- p 2 , v ((3, t, s -2/("+D ), 



-(i/+l)/2 
2,1/ j 



i<-<^r i)/a , 



-i. 



s > 6 



-(I/+1J/2 

1,1/ J 



s>0. 



Remark 2.5. The spectral measure corresponding to (2.9) is concen- 
trated on [b 2 [T + 2 j &i [r + 2 ]> having a density on (6 2 [T ; &i ^ +1 ^ 2 ) and 
atoms at each of the boundary points of the interval. 



Proof of Proposition 2.3. 
-\ogP{Z(0)< Wl ,Z(t)<w 2 } 



irvT(v/2) 



< (1 + f3 2 u 2 /u)-^ +1 )l 2 (1 + (3 2 {u - t) 2 /u)-^ +1 )' 2 , 
max^ , } du. 

-oo 



U'2 



Now 



( 1+)3 2 M 2 /ly) -.(H-l)/2 (l +/ 32( M _ t )2/^)-(^+l)/2 



if and only if (1 — x)u 2 — 2tu > 



— ^ — i 2 . This is equivalent to 



(2.11 



1 -x 



> 



(1-x) 2 [3 2 



if < x < 1 and to the reversed inequality if x > 1. Hence if < x < then 
(2.11) holds for all it and we get the last line of (2.9). Similarly if x > b 2tU , we 
get the first line of (2.9). In the case b\ )U < x < 1 condition (2.11) becomes 



t 

1 — x 



> 



xt 2 V 
(1-x) 2 ~J 2 
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that is, 



u> 



l-x 



+ 



xt 2 



or u< 



xt 2 



(l-x) 2 f3 2 'l-x \(l-x) 2 13 

Hence the integral over this interval becomes 

PT((v + 1) j2) rt/(i-x)-y/xty(i-x)*-v//3* (x + pp u 2/ u y0>+i)/2 



2 ' 



■kvT{v/2) 
/3T((i/ + l)/2) 



du 



+ 



(l + ^V/i/)-^ 1 )/ 2 



iruT(v/2) Jt/(i-x)+yjxt 2 /{i-x) 2 -v/i3 2 
0t 



du 



10l 



+ 1 - ^,1 1 



1 — X 



xt 2 V 

(1-X)2 "^2 

+ /3 



xt 2 



(l-x) 2 /? 2 



This, in combination with the integral stemming from the case 
V+P« a M- lv+1) » < (i±femZ^ ; gives the fourth h ne of ( 2 . 9 ). Simi- 
larly in the case 1 < x < 62,1/ we get the second line of (2.9). Note that in the 

if and only if u < t/2 for 



(l+/? 2 ix 2 / I /)-("+ 1 >/ 2 > (l+f3 2 (u-t) 2 /v)-W>/ 2 



W-2 



case x = 1 , 

t > and u > t/2 for t < 0. Hence we get for (2.9) £-P{T Vjl < 0\t\/2}. □ 
Next we move to the two-dimensional models. 

Proposition 2.4. For t = (ti,t 2 ) G M 2 and the exponential model, for 
wi,w 2 > 0, 



logP{Z(0,0)<w u Z(t ll t 2 ) < w 2 } 



(2.12) 



1 



1 

1 1 

1 1 1 

1 



1 1 



-/3(|ti| + |t 2 |)/2 



-/3(|ti| + |i 2 |)/2 



-/3(|ti| + |t 2 |)/2 



(wi,w2) g Ai, 

4 4 V^2 

(101,102) G A 2 , 

1 + p uM\t 1 l\t^ 



W± W 2 y/WlW2 



(wi,w 2 ) G A 3 , 
1 + /3 N+JM + Imf^ 



V101 



1 

102 ! 



(101,102) G A 4 , 
(101,102) G A 5 , 
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where 
Ai = 

A 2 = 

A 3 = 

A 4 = 

A 5 = 



P 



U'2 



(wi,w 2 ) : -In — <-(|*i| + |t 2 |) , 



(wi,w 2 ) 



(\h\ + M) < -in (^j < \h\ A \t 2 \ - 1*1 1 V |t 2 ||, 

<|<l|V|t 2 |-|t 1 |A|t 2 |k 



1 



K,^ 2 ):|tl|A|t 2 |-|ti|V|t 2 |<-ln 

(ioi,ti72) : |ti| V |t 2 | — |*i| A |t 2 | < -sin 

P 



Wi 

w 2 
w 2 



<N + |*2| 



(wx,w 2 ) 



p 



111 



w 2 



Proof. We work out the integral as in the one-dimensional case on 
the areas of u defined by \ui — ti\ + \u 2 — t 2 \ — \u\\ — \u 2 \ ^ ^ln(^). It is 
convenient, similarly to the one-dimensional case, to consider separately the 
nine areas defined by the position of u\ with respect to and t\ and by the 
position of u 2 with respect to and t 2 . So, for example, if < t\ < t 2 , we 
can write 



l^i - h\ + \u 2 - t 2 \ 



\ u i\ 



\ u 2\ 



—tl — t 2 , 

—2u 2 + t 2 —tx, 
t 2 — h, 

—2ui + t\— t 2 , 
-2ui - 2u 2 + ti + t 2 , 
-2ui + ti + t 2 , 
*i —t 2 , 

-2u 2 + t 2 + ti, 
[h+t 2 , 



if u± > ti and u 2 > t 2 , 
if ui > ti and < u 2 <t 2 , 
if ui > ti and u 2 < 0, 
if < u\ < ti and u 2 > t 2 , 
if < u\ < ti and < u 2 < t 2 , 
if < u\ < ti and u 2 < 0, 
if u\ < and u 2 >t 2 , 
if ui < and < u 2 < t 2 , 
if u\ < and u 2 < 0. 



The calculations are complicated but not difficult. □ 

Proposition 2.5. For t = (t\,t 2 ) £ R 2 and the normal model, we have 
for wi, w 2 > 0, 



(2.13) 



\o g P{Z(0,0)<wi,Z(ti,t 2 )<w 2 } 



Wl 



|t|/? 



+ 



1 



w 2 



t\p 



Wl 



1 



log — H $ 



w 2 



IMP, 1 1 w i \ 

2 |t|/3 w 2/ ' 



i/iat is, it is the same as in the one- dimensional case with \t\ replaced by |t|. 
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— /3^|u|^/2 — 1|"^/2 

Proof. As in the one- dimensional case. Now > if 

Wi — W2 

and only if f3 hu ^ 2U2 < ^ + ^ log ^. Next note that if the vector (Ui,U 2 ) 
has a standard two-dimensional normal distribution, ^i+^a is also stan- 
dard normal. The rest of the proof is as in the one-dimensional case. □ 



Proposition 2.6. Fort = (h,t 2 ) £l 2 and the t-model, for w\,w 2 >0, 

-logP{Z(0,0) <wi,Z(h,t 2 ) <w 2 } 
(I 

— , 0<w 2 <b 2 ^w 1 , 
w 2 

— P{{Ti,T 2 ) E A 1>a } + —P{(T U T 2 ) € A^J, 
W\ w 2 

b 2y a w l <W 2 <W\, 

(2.14) = < -P{T X < |t|/2}, Wl = w 2 =:w, 

w 

— P{(Ti,T 2 ) £A c la } + -Lp{(Tx,T 2 ) G A 2 , a }, 
W\ w 2 

w 1 <w 2 < b^wi, 

where (Ti,T 2 ) is a random vector with bivariate t-density (1.7), 
/3 2 |t| 2 



i + 



b 2 , a = i + ] — + 



4(a-l) v/2(a-l) 
/?|t| 



21 + 12 



'1 + 



/? 2 |t 



21 + 12 



4(a-l) ^2(a-l) 



'1 + 



/? 2 |t 



(a-1)' 



.4 



(«l,« 2 ) GM 2 : Ui 



tl 



1 — x 



+ li2 



*2 



1 — X 



M, a = { (ui,u 2 )GK 2 : (tti 



< 

£ix 
1 -x 

< 



(1-x) 2 



x|t| 2 2(a-l) 

1 — X 

2(a- 1) 



+ I «2 

xltl 2 



(l-x)s 



and 



x ■ 



w 2 ) 



1/a 
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Proof. Analogous to the one-dimensional case. □ 

Proposition 2.7. For t = (ti,t 2 ) G M 2 and the general normal model, 
for wi, W2> 0, 

-logP{Z(0,0) < Ti>i,Z(ti,t 2 ) < w 2 ] 



(2.15) 



where 



1 



1 / Vt T E-H 



W 2 \ 

: log — I 



1 /Vt^E-H 

+ — * + 

t»2 \ 2 



Vt T E~H w 2 )' 



s- 1 



0? -p&ft 
-PP1P2 Pi 



Proof. As in the one-dimensional case. Now 



-u T £ _1 u/2 



> 



-(u-t) T £ -1 (u-t)/2 
U)2 



if and only if u' J 'E~ 1 t < - +log^. Next note that 
if the vector U = (C/i, C/2) has a bivariate normal distribution with mean 
value zero and covariance matrix E, U T E _1 t has a normal distribution 
with mean value zero and variance t T T,~ l t. The rest of the proof is as in 
the one-dimensional case. □ 

3. Estimating the dependence parameter (3. We consider a sequence 
of independent, identically distributed stochastic processes with continuous 
paths 

{Xi(t)} tm , i = l,2,.... 

We assume that the processes are in the max-domain of attraction [as pro- 
cesses in C(R)] of a max-stable stationary process {Z(t)} t £m. such that the 
related process Z (see the Introduction) has exponential spectral function 
(2.2) or (2.12), or normal spectral function (2.6) or (2.13), or t spectral func- 
tion (2.9) or (2.14) as discussed in Section 2. For definition of convergence 
and convergence criteria see [4]. 

Ideally it would be nice if we could assume that we have observed the 
sample paths of n processes X as a basis for estimation of the main parame- 
ter j3. However, in reality this is too much to expect. Usually one can observe 
the n processes only at finitely many points in space, say t\,t2, ■ ■ ■ , id- 

In this setup we propose estimators for (3 that are closely related to an 
extension of the estimator R(x, y) for the dependence function R(x, y) which 
was introduced by Mason and Huang (see Huang [9]), 

1 n 

R tl ,...,t d ( x ^ ■■■,x d )-= £ Xl / {-Y I (ti)>x n _ [fea;i]+li „(i 1 ),...,x l (i (i )>x n _ [fe , d]+li „(t d )} 
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with k -C n, where {Xi tn (tj)}f =1 are the nth-order statistics of {Xi(tj)}f =1 
for j = l,2,...,d. 

Let us now introduce our estimators in R. All three densities have the 
form (3ipo([3t). Corollary 3.2 below states that 



i^i td (l,...,l) = 2 / Mt)dt 

J (f3/2)(max 1 <j<dtj—mm 1 <j<dtj) 
2{ 1 — F ( — ( max tj — min t,-, 

\2 \l<j<d J l<j<d J 



with F(t) := f_ 00 ipo(s)ds. It follows that 





2F-(l-i? tl>t2 ,... jtd (l ) l,...,l)/2) 



max i<i<d*i ~ mm i<j<d*i 
Hence we introduce the estimators 

^(i) . = 2F-(l-^ 1 , f2 ,,,,, td (l,l,...,l)/2) 
max i<j<<i*j - mini<j<d*j 

and 

A := ith) K £ <d fe^u F ^ (1 " 1)/2> ' 

Considering the two-dimensional space, note that the standard normal and 
Student densities are spherical symmetric. This means (cf. proof of Proposi- 
tion 2.5) that it is sufficient to consider the marginal distribution and hence 
we introduce the estimators 

with F the (common) marginal distribution corresponding to (fo (standard 
normal or Student). The exponential model in two-dimensional space is more 
complicated. We consider 

2 

' 2 := did- Y) ^ ( ^ a i,^fio, m { R tj,tm( 1 A)) 

with aj jTn := — t£ |, the absolute difference of the first components of 
tj and t m , and 6 Jim := |tjj ~~ ^2 I an< ^ the i nverse function of ^(1 + 
^ min(a, b)) exp{— ^(a + 6)}, which is decreasing in (3 for a, 6 > 0. 

Note that f3^ is simpler than (3 and summarizes the information of the 
sample in a somewhat more crude way. We could not find analogues of /3« 
in two-dimensional space since we were unable to calculate explicitly the 
necessary higher-dimensional distributions. 
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All the mentioned estimators are consistent and asymptotically normal 
under the appropriate conditions. We now state the results. The proofs will 
be given at the end of the section after some lemmas. First we consider 
consistency. 

Theorem 3.1. Suppose that the normalized sequence of maxima [see 
(1.1)] converges weakly to {Z(t)} in C{ — oo,oo). For sequences k = k(n) — > 
oo, k(n)/n — > as n — > oo (recall that the sequence k figures in the definition 
of R), all the indicated estimators are weakly consistent for [3. 

Also the estimators are asymptotically normal under some extra condi- 
tions, that is, Vk(j3 — (3) has asymptotically a normal mean zero distribu- 
tion. In order to describe the asymptotic distribution more accurately, we 
now state a slight extension of a result of Huang and Mason. 

Proposition 3.2 ([9], pages 29 and 43). Let {(X^h), . . . ^X^td)}^ be 
i.i.d. random vectors with distribution function F. Suppose that the marginal 
distributions F{ (i = 1,2, . . . ,d) are continuous and strictly increasing. De- 
fine 

F tlr ..,t d (xi,- ..,x d ):=l- F(Fr(l -x 1 ),...,FT(l-x d )) 

= P{1 - F 1 (X(t 1 )) < Xl or ■■■ or I- F d (X(t d )) < x d }, 

where the arrow denotes inverse function. Suppose that for allx±, x%, . . . ,x d > 
0, Xi + X2 + • • • + x d > and a positive function L, 

(3.1) lim^F tut2 ,...,t d (txi,tx 2 , ■ ..,tx d ) = L tu t 2 ,...,t d {xi,x 2 , . ..,x d ). 

Next we introduce the definition of a function R which is connected with 
the function L as follows. Let v^fa t d be the measure that satisfies for 
x 1 ,x 2 ,...,x d >0 

vti,...,t d {( s i'---i s d)\ s i < x i or ••• or Sd<Xd} = L tl ,...,t d (xi,...,x d ). 
Such a measure exists by virtue of (3.1). The function R is given by 

Rt 1 ,t 2 ,--;t d ( x i,X2,---,x d ) ■■=vt 1 ,t 2 ,...,t d ([0,xi] x [0,x 2 ] x ••• x [0,x d ]). 
Define 

1 n 

Rt lt ...,t d ( X l' ■■■> X d)-= ^X]' f {^i(*l)>^r,-[* a : 1 ]+l,n(*l)."- ) Jfi(*d)>-y»-[fcx ( ,]+l,n(*d)}' 

i=l 

where X n _[ kx .] + i n (tj) is the (n— [kxi\ + l)st-order statistic of X\(tj), X 2 (tj), 

. . . , X n (tj), j = 1,2, . . . ,d. Then for all x±, x 2 , . . . , x d > 0, x\ +x 2 -\ \-x d > 

0, 

(3.2) R tl ,t2,...,t d (xi,x 2 ,...,x d ) ->■ R tu t2,...,t d (x 1 ,x 2 ,...,x d ) in probability, 
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n — > oo, k = k(n) — > oo, k/n — > 0. 

Suppose now that the function L has continuous first-order partial deriva- 

(i) 

tives Ll l t2 t , j = 1, 2, . . . ,d. Moreover suppose that for some a > 
F tl ,t 2 ,...,t d {tx 1 ,tx 2 , ■ ■ ■ ,tx d ) 

(3.3) 

= t{L tl>t2 ^ td (xi,X2, ...,x d ) + 0(t a )}, 1 1 0, 

uniformly on xf + x^ \-x 2 d = 1, x% > 0, % = 1, 2, . . . ,d. Then for a sequence 

k = k(n) — ► oo with k = o{n 2a ^ 2a+l ^), n — > oo, 

(3.4) \fk{R tl ^ td (xi, . ..,x d )- Rt u ...,t d ( x y> ■ ■ -> x d)} B tl! _ ttd (xi, . ..,x d ) 
in D(WV), where 

B tl ,t 2 ,...,t d (xi,x 2 ,. . . ,x d ) 

= W tljt2r .. M (x 1 ,X2,...,X d ) 

- L tl,t 2 ,...,t d ( x i> x 2, x d )W tl , t2 ,...,t d (xi,0, . . . , 0) 

L ti,t 2 ,...,t d ( x i, x 2, • • • > x d)W tl ,t 2 ,...,t d (0, 0,...,x d ) 

and W is a continuous mean zero Gaussian process with covariance structure 

EWt 1 ,...,t i (x ll x d )W tl ,..., td (yi, ...,y d ) = v tl ,...,tMxu-,xd n A vu-m) 
with 

A XuX2t ..., Xd :={(t 1 ,t 2 ,...,t d )\t 1 <x 1 ort 2 <x 2 or ■■■ ort d <x d }. 

The function L is called the tail dependence function and it is directly 
related with the extreme-value limit distribution. In fact if F is in the do- 
main of attraction of an extreme- value distribution G, condition (3.1) holds 
with L(xi,...,x d ) = -log G((- log Gi)^(a;i ),..., (-log G d )*~(x d )) where 
Gx, . . . , G d are the marginals of G. 

dimensional space 
.,x d e-^ s - td \)ds 

.,x d e-^ s - td \)ds. 



Remark 3.1. For the exponential model in one- 



-P\s~t! 



Lt u ...,t d {xi,...,x d ):=- / max(xie 

^ J — oo 

and (see Lemma 3.1 below) 

P f°° 

Rt u ...,t d {xi,...,x d ):=- / min^ie - ^ 8- * 1 ', 

J— oo 

Similarly for the normal model and the i-model. 
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Theorem 3.3 (One-dimensional space). Suppose for the processes 
{X n (t)} the limit relation (1.1) holds and moreover the second-order condi- 
tion (3.3): for some a > 

F tl> ,.. jtd (t Xl , . . .,tx d ) = t{- logP{Z(h) < 1/xi, . . • , Z{t d ) < l/x d } + 0(t a )}, 

1 10, uniformly on x\ + x\ H V x d = 1, Xi > 0, i = 1, 2, . . . , d. 

Ifk = k(n)^oo, k(n) = o(n 2a ^ 2a+l ^) , asn^oo, then 

(3.5) Vk~0i - /?) - ^TJ E ^ - " * m |/2) 

and 

^fc(/3 (1) -/3) 
(3.6) 

2^ 1 , t2 ,„„ fd (l,l,...,l) 1 

maxi<j<d*j - ™dni<i<d*j <£>o03(maxi<.,< d tj - min 1 < J < d t J )/2) 
in distribution. Here B is as in Proposition 3.2. 

Theorem 3.4 (Two-dimensional space). Under the same conditions, 

i- * 2 „2B t t (lA) 1 

(3.7) ^(ft-fl^ £ 



and 

(3.8) ^0 e ,2 " /3) - E ^,t m (l, Wa jm!bjm (Qa jm , bjm W) 

^ ' j<m 

in distribution, where a jm := — ti |, 6j m := |#2 — ^2™^ I- 

For the estimation of the general normal model we proceed as follows. 
Write 



ST*- 1 



1-P 2 



/?i 2 



2 

2 



with — 1 < p < 1 and /?2 > 0. For two sites tj and t. m (1 < j < m < d) we 
write 

(tj t m ) S (tj t m ) 
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where 



and 
(3.9) 

Now note that 



t j = (tf\tf ) f for j = 1,2,... ,d, 

{tf-t%f 
(tf-t^)itf-t^) 
(tf-t^f 



PI 



a := 



1-p 2 



2 - i? tj , tm (l, 1) = L tj>tm (l, 1) = 2$( v /|t J - - t m rs-i|t, - t m |/2). 
Hence we define estimators 

(3-10) 4, m := (2*^(1 - R tj ,t m (l, l)/2)) 2 . 

Using the result of Proposition 3.2 and Cramer's delta method we get 



(3.11) y/k(Qj, m - tj >m a) 4 2B tj>tm (l, l)y / t^(«/»( v /t^/2))- 1 , 

where <f> is the standard normal density. Now compose the (d(d — l)/2)- 
dimensional vectors 

Ql,2 
Ql,3 

q: = 



-Q 



and 



d-l,d- 



J l,2 



J l,3 



Then 
(3.12) 



Vk(q - Ta) 4. b 
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with 



s tl)t3 (ia) A /^(0( A A!>/2))- 1 



^t d _ ljtd (l, l)V / C 1 , d a(^( v /tJ_ 1(i a/2))" 1 
Next define 

(3.13) a:= (^r)" 1 ^^. 
Then 

(3.14) Vk{k - a) = (r T r)- 1 r T >/fc(q - ra) S (r^r)- 1 ^. 

By solving the equations (3.9) we get 



and 



Pi- 

ft: 

P 



1 0,1 



4a 3 ' 



'a 3 



4ai 



a 2 



2^aia 3 



Cramer's delta method now gives the joint asymptotic normality of the 
estimators 



(3.15) 

(3.16) 
and 
(3.17) 



Pi: 
fa' 

P 



'Oi 



«2 



4a 3 ' 



'a 3 



;,2 



4ai 



a 2 



2v / «ia3 

For the proofs of the theorems we need a number of auxiliary results. 

Lemma 3.1. Suppose for some measure v on R d and some positive in- 
tegrate functions g\ , . . . , g d 



u([xi,oo) U • • • U [x d ,oo)) 



max{5i(x - xi), . . . ,g d (x - x d )} dx. 
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Then 

/oo 
mm{g 1 (x - xi), . . .,gd(x - x d )}dx. 
-oo 

PROOF. Note that 

u([xi,oo) U [x 2 , oo) U • • • U [x d ,oo)) 
d 

= u (l x i' °°)) + H v ([ x h°°) n i x ji °°)) 
i=l i^j 

h h (-i) d_1 i/([^i,oo) n [x 2 , oo) n ■ • • n [x d , oo)), 

and for any real a\, a 2 , ■ ■ ■ , a d , 

d 

max(ai, a 2 , ■ ■ ■ , a d ) = Y~] a? + (—1) Y~] o-i A eij + • • • + (— 1) _1 ai Aa 2 A - ■ ■ Aa d 



i=l 



(both follow by induction). Then the statement also follows by induction. 
□ 

Corollary 3.1. Hence for our models in R 
Rt 1 ,t 2 ,...,t d (xi,x 2 , ...,x d )-- 



mm as. 

oo 1 Xi 



Remark 3.2. We apply Lemma 3.1, for example, to the functions gAx) :- 

e -P\A 



Lemma 3.2. Suppose p is a probability density on R, p(x) =p(—x) for 
x > and p(x) is decreasing for x > 0. Then 



(3.18) 



min{p(|s - ti|), . . . - t d \)} ds 



(l/2)(maxi<j<d%— mini<j< d %) 



p(s) ds. 



Proof. Note that 
min{p(|s - ti|), . . . ,p(\s - t d \)} 



mm 



K 



min 



l<j<d 

s — min t 



■P 



max t, 
i<j<d J 



l<j<d 

max f 

l<j<d 



s > min tj + k ( max i 



l<j<d 



min i 



i<j<d i<?<<2 



s < min tj + 2 I max t 



min i 



i<j<d i<i<<i 
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So the integral on the left-hand side equals 



P 

mini< J < d f J + (l/2)(max 1 < j < d tj-min 1 <j< d t j ) 



S — min tj 

l<j<d 



ds. 



□ 



Corollary 3.2. Hence for our models in R 

/■oo 

R tu ... ttd (l,...,l) = 2 0(3) ds 

J (1/2) (max, < 1<0 rt, -mini < ,< rf t, ) 



'(l/2)(max 1 < j < d tj-mini< J < d t :l ) 

= 2/ 1 — -F^— f max tj — min 

I V 2 \i<j<d J i<j<d J 

Lemma 3.3. Under the conditions of Theorem 3.3, 

Vk[R tl ,..., td (h ■ ■ ■ , 1) - 2 [l - F(| (max - mm t,-) ) 
(3.19) 

->• -Bti,...,t d (l, ••-,!) 
in distribution with B as in Proposition 3.2. 

Proof. Combine Proposition 3.2 and Corollary 3.2. □ 

Lemma 3.4. Under the conditions of Theorem 3.4, for the standard nor- 
mal and Student models, 



1 



(3.20) \/fc{i? t ., tro (l, 1) - 2(1 -F(£(\tj - tm|))) } - £ t ,, tro (l, 
m distribution. For the exponential model 

^{«t 3 ,t m (l, l)-|(l + f min(|t? - 4 m) |,|4 i} - 4 m) I) 

(3.21) xe -(/V2)(i^M m) i+i4M m) i) 

-S t ., tra (l,l) 
in distribution, where tj = (t^ ' , t^), t m = , 4 )• 

Proof. Propositions 2.5, 2.6 and 2.4 give Lt, t m (l,l), which is 
2-^.^(1,1). □ 

Proof of Theorem 3.1. Follows immediately from statement (3.2) of 
Proposition 3.2 and Lemma 3.2. □ 

Proof of Theorems 3.3 and 3.4. Immediately from Lemmas 3.3, 3.4 
and Cramer's delta method. □ 



SPATIAL EXTREMES 



23 



Acknowledgments. Comments of two referees made us aware of Smith 
[14] and helped greatly to improve the presentation. 



[1] Balkema, A. A. and de Haan, L. (1988). Almost sure continuity of stable mov- 
ing average processes with index less than one. Ann. Probab. 16 333-343. 
MR0920275 

[2] Brown, B. M. and Resnick, S. I. (1977). Extreme values of independent stochastic 
processes. J. Appl. Probab. 14 732-739. MR0517438 

[3] Coles, S. G. (1993). Regional modelling of extreme storms via max-stable processes. 
J. Roy. Statist. Soc. Ser. B 55 797-816. MR1229882 

[4] de Haan, L. and Lin, T. (2001). On convergence toward an extreme value distribu- 
tion in C[0, 1]. Ann. Probab. 29 467-483. MR1825160 

[5] de Haan, L. and Pickands, J., Ill (1986). Stationary min-stable stochastic pro- 
cesses. Probab. Theory Related Fields 72 477-492. MR0847381 

[6] de Haan, L. and Resnick, S. I. (1977). Limit theory for multivariate sample ex- 
tremes. Z. Wahrsch. Verw. Gebiete 40 317-337. MR0478290 

[7] Eddy, W. F. (1980). The distribution of the convex hull of a Gaussian sample. J. 
Appl. Probab. 17 686-695. MR0580028 

[8] Falk, M., Husler, J. and Reiss, R.-D. (1994). Laws of Small Numbers: Extremes 
and Rare Events. Birkhauser, Basel. MR1296464 

[9] Huang, X. (1992). Statistics of bivariate extreme values. Ph.D. dissertation, Tinber- 
gen Institute. 

[10] Husler, J. and Reiss, R.-D. (1989). Maxima of normal random vectors: Be- 
tween independence and complete dependence. Statist. Probab. Lett. 7 283-286. 
MR0980699 

[11] Pickands, J., Ill (1981). Multivariate extreme value distributions (with discussion). 

Bull. Inst. Internat. Statist. 49 859-878, 894-902. MR0820979 
[12] Schlather, M. (2002). Models for stationary max-stable random fields. Extremes 5 

33-44. MR1947786 

[13] Sibuya, M. (1960). Bivariate extreme statistics. I. Ann. Inst. Statist. Math. 11 195- 
210. MR0115241 

[14] Smith, R. (1990). Max-stable processes and spatial extremes. Unpublished 



REFERENCES 



manuscript. 



Econometric Institute 

Erasmus University 

Rotterdam 

Netherlands 

E-MAIL: ldehaan@few.our.nl 



Department of Statistics 



and Operational Research 
Faculty of Science 
University of Lisbon 
Campo Grande, Bloco C6, Piso 4 
1749-016 Lisboa 
Portugal 

E-MAIL: tpereira@fc.ul.pt 



